Dynamics of pair correlations in the attractive Lieb-Liniger gas 
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We investigate the dynamics of a ID Bose gas after a quench from the Tonks-Girardeau regime to the regime 
of strong attractive interactions applying analytical techniques and exact numerical simulations. After the 
quench the system is found to be predominantly in an excited gas-like state, the so-called super- Tonks gas, 
however with a small coherent admixture of two-particle bound states. Despite its small amplitude, the latter 
component leads to a rather pronounced oscillation of the local density-density correlation with a frequency 
corresponding to the binding energy of the pair, making two-particle bound states observable in an experiment. 
Contributions from bound states with larger particle numbers are found to be negligible. 

PACS numbers: 03.75.Kk, 67.85.De, 05.30.Jp. 34.20.-b 



Ultra-cold quantum gases in reduced spatial dimensions 
have attracted a lot of attention in recent years This is 
because on one hand quantum effects play an increasing role 
in lower dimensions and on the other hand these systems be- 
came experimentally accessible using ultra-cold atomic gases. 
A striking example for the pronounced effects is the effective 
fermionization of a one-dimensional Bose gas with repulsive 
interactions, described by the Lieb Liniger (LL) model 101, 
leading to the so-called Tonks-Girardeau (TG) gas Isl-IHt]. In 
the attractive case the ground state of an iV-particle LL gas 
is the highly localized McGuire cluster state |^. In the ther- 
modynamic limit the gas is unstable, preventing direct exper- 
imental studies of ground-state properties. However, in dy- 
namical setups, attractive and repulsive gases are equally well 
accessible. A recent milestone in this direction is the creation 
of the super Tonks-Girardeau (sTG) gas ||7[lll] by Haller et al. 
Jgtl realized by a rapid sweep through a confinement induced 
resonance from the strongly repulsive to the strongly attrac- 
tive side. The sTG gas is a highly excited, gas-like eigenstate, 
that does not contain any bound pairs or higher particle clus- 
ters. We here analyze the dynamics of this quench process 
by numerical simulations employing the time evolving block 
decimation (TEBD) ifloilTlll algorithm recently applied to the 
simulation of the relaxation dynamics of the repulsive LL gas 



1 12[] and a number of lattice models Ill3l4l5 1 



The Hamiltonian of a one-dimensional, trapped Bose gas 
with local interactions is given by the Lieb-Liniger model yfl 
with an additional potential term 

(1) 

in units were m = 1. gis, the interaction strength, which 
is related to the ID scattering length via g = — 2/aiD, 
and can have both signs. It is often characterized by the di- 
mensionless Tonks parameter 7 = g/g. V{x) = ^w^x^ de- 
scribes an harmonic trap. 

The spectrum of ^ is quite different depending on the sign 
of the interaction, however the positive- and negative-7 spec- 
tra agree in the limits of weak as well as strong interactions 
1 1611 . This can be most easily understood from the homoge- 



neous problem of two particles with periodic boundary condi- 
tions (PBC): Fig.[T]shows the lowest lying states with vanish- 
ing center-of-mass momentum. The non interacting ground 
state (left and right end of the figure at energy = 0) has a 
constant relative wave function between the two particles. For 
finite interactions, the wave function must obey contact con- 
ditions, equivalent to the interaction Hamiltonian \2^. It de- 
velops a peak at zero inter particle distance when an attractive 
interaction (7 < 0) is turned on and eventually forms a closely 
bound pair with binding energy 7^ (see below). On the repul- 
sive side (7 > 0) a dip-like kink emerges with increasing in- 
teraction, which eventually makes the wave function vanish at 
coinciding particle positions - this is the famous fermionized 
TG gas. Most importantly when approaching the strong in- 
teraction regime from the attractive side, the first excited state 
adiabatically connects to exactly the same fermionized state - 
the sTG gas. This matching continues for higher excited states 
and the scheme can be generalized to many particles, where 
the bound states can be classified by the number of dimers, 
trimers, etc. 




FIG. 1 : (Color online) Spectrum of LL Hamiltonian for two particles 
with vanishing total momentum on a ring of length i = 1 as a func- 
tion of inverse interaction strength. One recognizes equivalence of 
the spectra at vanishing (I7I — 0) as well as infinitely strong inter- 
actions (I7I — )■ 00). 

Fig. [T] indicates, that a quench from the TG regime to the 
strongly attractive regime will put the gas to good approxima- 
tion in the sTG state. Experiments have successfully demon- 
strated this, while the difference between TG and sTG can be 
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detected by their different compressibility 

We here consider a gas of N particles confined by a har- 
monic trap initially being in the ground state for 7 = +00. At 
t = the interactions are switched to the strongly attractive 
side 7^—1. The trap plays a minor role, since interactions 
give the relevant time scale dominating over the kinetic energy 
resulting from the trap confinement. 

We simulate the full many-body dynamics using an exact 
numerical technique. To this end the continuous model ^ 
is discretized, resulting in the sparsely filled Bose-Hubbard 
model i fisil . The lattice is actually finite, but comprises all of 
the gas (which does not expand globally on the time scale in 
question). The dynamics is then simulated using the TEBD 
scheme, employing a fourth order trotter decomposition lfl9ll . 
For the specific setup it is necessary that the conservation 
of the total particle number is taken into account explicitly. 
While time dependent simulations are generally limited to 
short times due to the linear growth of entanglement entropy 
I 20I , 21 1, this it not crucial here. Although we also observe 
such a linear growth, the increase is slow since we are close 
to an eigenstate. Thus we can go much beyond the time scale 
of interactions as in (l^ for as much as iV = 18 particles on 
a 1280 sites lattice using a rather small matrix dimension of 
X = 100 in the algorithm. 

Fig.|2]shows the dynamical evolution of the density-density 
correlations, where we fix one position at the center of the 
cloud. Time is given in all figures in units of A/g^. The 
initial state shows the typical feature of fermionization, i.e. 
f^^^-* (0, 0) is zero and rises to one (no correlation) on a length 
scale proportional to the average inter-particle distance. In the 
limit 7 — 00 the correlations do not show much resolvable 
dynamics on the scale of the figure because the initial TG is 
close to the sTG state. However for the moderate interaction 
strength chosen, we see g^^-' rising sharply around zero dis- 
tance. This must be due to transitions to states other than the 
sTG state. One finds that the characteristic length scale of the 
peak at the origin is given by the ID scattering length aid. 
This indicates a finite admixture of the N = 2 McGuire clus- 
ter state. Note, that since / dx{i/'^ {x)i''^ {y)ii {y)ii (x)) = 
N — 1 the integral over g*^^) (0, x) must be constant in time, as 
long as the density can be assumed to be homogeneous. The 
increase at xi — 2:2 =0 must therefore be accompanied by 
a decrease at larger distances as seen in the correlation waves 
building up in Fig.|2] The insert of Fig. |2]shows another inter- 
esting feature: Apart from small distances, where oscillations 
continue, the correlations become quickly stationary and show 
rather good agreement with g^'^^ (0, x) of a hard sphere Tonks 
gas (HS) gas with hard-sphere radius aiD- This verifies a re- 
cent finding by Girardeau and Astrakharchik ll22ll . who have 
shown that the wave function of the sTG gas is identical to that 
of a TG gas where a small avoided volume of size uq « oid is 
inserted into the two-particle wave function around zero dis- 
tance (HS Tonks gas). In Fig|2l the density-density correlation 
of the HS Tonks gas is approximated by g^^^ of a TG gas on 
a reduced volume where the hard-sphere volume is excluded. 
This corresponds to the t — curve multiplied by 1 — aiD£' 



and shifted by aiD. For separations larger than the interparti- 
cle distance 1/ g the non local nature of the TG - HS mapping 



*Hs(a:^i < X2 < ■■ ■ < xn) 
= ^-TG ixi,X2 +ao,...,XN + {N ~ l)ao) 

prevents however a simple calculation of g'^' (0, x). 



(2) 




FIG. 2: (Color online) Time evolution of the non-local density- 
density correlation function for interaction strength 7 — —18.7931 
calculated for A'' = 18 particles using TEBD. Curves are shown for 
(black), 1/4 (dotted), and 1/2 (red, maximum value at 2; = 0) os- 
cillation periods of the local correlation, insert: g'^' at times of 3.5 
(black), 4 (dotted), and 4.5 (red) periods. The vertical dashed line 
indicates x = aio. The blue dashed line shows an approximation to 
the HS ground state mentioned in the text. 

Fig.[3t shows for equal positions as a function of time 
for various values of the interaction strength. The correlation 
function grows as a power law with exponent depending on 
the interaction strength. By a linear fit to the numerical data, 
we find it growing from 1 in the free case to a value of about 
4/3 in the strongly attractive case. We see that g*^^' rises up 
to a finite value much smaller than 1 for reasonably strong in- 
teraction. This reflects the fact, that most of the gas ends up 
in the fermionized sTG state. For longer times and stronger 
interactions, we observe however a rather peculiar oscillatory 
behavior with large modulation depth. This can only be un- 
derstood as a result of an interference between two compo- 
nents of a coherent superposition, the sTG state and a bound 
state. In fact the oscillation frequency coincides with the bind- 
ing energy of a pair of particles in the McGuire state. Thus 
the dynamics is strongly affected by the contribution of bound 
pairs. Moreover, there is no sign of a relaxation, as observed 
in the repulsive case [l2l. In Fig. [3}' the local three-particle 
correlation g^^^ is plotted. One recognizes that g'^^'> remains 
extremely small, showing that higher-order cluster states are 
not formed in the interaction quench. This agrees well with 
the finding in 12211 where the overlap of the TG wave function 
with the McGuire cluster state was calculated. 

We will show now that the two-particle correlations in the 
gas can be very well reproduced by a system containing only 
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FIG. 3: (Color online) Time evolution of the local two particle den- 
sity density correlation in a system of A'^ = 18 particles calculated 
via TEBD (colored lines), gray line: two particle case with PBC 
and 7 = 0. insert: Time evolution of local three-particle correlation 
g'^' (red dashed) at the trap center and g'^' (black) for comparison 
at 7 = —145. (The artifacts for very short times are due to the finite- 
size time steps used by the numerical algorithm.) 



N = 2 particles. This is due to the fact that for strong interac- 
tions, eigenstates are well approximated by pair product states 
of the Jastrow-Bijl type 112211 
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We impose periodic boundary conditions (L = 1), which is 
reasonable for the comparison to the trapped gas, since it will 
be homogeneous in good approximation over the length scale 
of some inter-particle distances. The PBC problem gives ana- 
lytical expressions and allows to extract the scaling of certain 
quantities with 7 in the strongly interacting regime. This prob- 
lem has been solved for attractive interactions in the appendix 
of the original paper by Lieb and Liniger . We will use their 
solution in the following. 

The Hamiltonian for the two particle problem reads in first 
quantization H = {d\ + 91) + ~ 2:2). All eigen- 

states of the LL model can be constructed from coordinate 
Bethe ansatz lEsll . In the primary sector (0 < cci < 0:2 < L), 
the solution is ip{y = xi — X2) = 2Ae^^ cos(|(j/— i)). 
A is an interaction-dependent normalization constant, and 5 
is related to the scattering phase shift = — 2tan^^((fc2 — 
ki)/2-j) via 9 = S/2 — tt. Note that e'j is not a simple phase 
factor as S will be imaginary for the two-particle bound state. 

We will now calculate asymptotic expressions for the bound 
state ipb, where 7 — ?> —00, as well as TG and sTG states 
ip±, where 7 — > ±00. For the bound state ipb we need to 
find an imaginary solution of the Bethe equation 6/2j = 
l/tan((5/4). Any S that is not purely real, must be purely 
imaginary, as 7 is real. Substituting S ~ id we find in the 
strongly interacting limit S = (—27). With this we calcu- 
late the normalization of the wave function, yielding ^■ 

iJs/2. Thus the wave function at coinciding particle posi- 
tions reads 



(^b(0,0) = 2Abe 3 cosh 
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(4) 



S/4 



As the density gis2 everywhere, this results in 

-7/2. 

We denote the lowest lying gas like states Lp± for 7 — > ±00. 
From the Bethe equation we see that a real solution S will 
be close to 27r. Expanding the tangent around its singularity 

at 7r/2 we get S ~ 27r ^1 — . For the normalization this 

means A± 1/ \/2, such that the wave function at coincid- 
ing particle positions becomes 



ip±iO,0) = 2A±e' 



u 7— — 00 
COS - > I 
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V2^. 

7 

(2) 
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The local two-particle correlation is in this case = 



The 7^^ scaling is well known 12411 and is the same as in the 
many particle case. 

The overlap between the initial TG gas state (/?o ~ 
lim-y^Qo and the bound state for finite 7 can easily be 
calculated. In the strongly interacting regime one finds e = 
(</Jol^fc> ^ -2x/2^7-3/2. 




FIG. 4: (Color online) Comparison of the many body results with 
the results from the two particle system with PBC in the case 7 = 
—89.0355. The thick black line corresponds to the two particle case, 
the thin green one shows the N = % particle case (calculated via 
TEBD), and the blue dashed one is the beating approximation l|9]l. 

We now want to calculate the local correlation dynamics in 
the two particle case. Numerically this is easily doable by 
constructing all eigenstates solving the Bethe equation and 
re-expressing the initial state using this eigenbasis. On the 
other hand from the above calculations we can derive simple 
approximations for the strongly interacting regime which are 
very good and give insight into the nature of the observed os- 
cillations. We can decompose the initial state |(po) according 
to 



(6) 



Note that for large I7I, is approximately normalized. One 
finds that for t = the wave function at coinciding particle 
positions is given by 



<^o(0,0,t = 0) 



% ^2^/2-. 

7 



(7) 



(po is however not an eigenstate, but is composed out of low 
lying gas like states, that have an energy spread much smaller 
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than the binding energy of the pair state. For small times one 
can ignore the energy differences and thus ignore the time de- 
pendence of (fQ. This results in the correlation 



r 



1 — cos(7^t) 



(8) 



This expression describes the initial increase of g'^^ as ob- 
served in the many-particle calculation (and in the exact so- 
lution of the N = 2 case) very well. It does predict, how- 
ever, oscillations with unity modulation depth, which is not 
true for the exact solution. The reason for this is that ipo con- 
tains in addition to the dominant, lowest gas-like state (i.e. the 
sTG state) also small admixtures of higher lying gas states 
which oscillate in time all with slightly different frequencies. 
For larger times these oscillations lead to an effective dephas- 
ing in the interference part of g'^' . On the other hand the direct 
contribution of these higher excited gas states to g^^^ is neg- 
ligible. Thus an approximation which is much better suited 
to describe the large time behavior is \ipo) « e\(fb) + If-)- 
Comparing Q and (|7]l shows, that we have only changed a 
factor of 2 such that 



{5-4cos[(7^ + 7r2)t]}^, 



(9) 



where we used that the sTG gas energy is tt^ for strong inter- 
action, giving a minor correction to the frequency. On short 
timescales this expression is of course invalid. However at 
times t > 1/7^ it becomes a much better approximation than 
(O as shown in Fig.|4] 

In summary we have shown by numerical TEBD simula- 
tions that an interaction quench of a ID Bose gas from strong 
repulsive to strong attractive interactions puts the gas predom- 
inantly into the lowest gas-like excited state, the sTG gas. 
There is however a small coherent admixture of two-particle 
bound states that results in a large amplitude oscillation of the 
local density-density correlation with a frequency correspond- 
ing to the energy difference between sTG gas and bound pair 
state. At the same time higher-order correlations remain ex- 
tremely small showing that more deeply bound, multi-particle 
cluster states are not formed in the quench. Analytical calcula- 
tions of the N ~ 2 case where shown to reproduce the results 
of the many-particle simulations with high accuracy. This in- 
dicates that the many -body state can be well approximated by 
a Jastrow-Bijl type pair product wave function, where each 
term is a coherent superposition of a gas-like state with a very 
small component of a two-particle bound state. The peculiar 
oscillations of g''^^ show furthermore that despite their small 
weight, the two-particle cluster states are accessible to exper- 
imental probes. For strong interactions, the bound pairs are 
highly co-localized. Since in all physical realizations of the 
LL model, the true inter particle potential is of finite range, 
details of the potential will show up in the binding energy. In 
this way, the two particle correlation dynamics can be used 
to measure details of the underlying true potential. In order 
to assess whether the effect predicted is accessible in current 



experiments let us estimate the beat frequency in SI units. 
Asymptotically it is given by a; = j'^hg^ /Am with m being 
the mass of the particle. This translates into uj = 7^-/Vw||/4, 
where is the longitudinal trap frequency, i.e., for the ex- 
periment presented in |2| the beating frequency is of the order 
of several kHz. This is comparable to the transversal trap fre- 
quency and thus an experimental observation requires either 
to go to smaller values of 7 or to use tighter traps. Finally 
our paper also shows that the TEBD algorithm is suitable for 
the simulation of dynamical processes in strongly interacting, 
continuous quantum gases. 
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